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Abstract. 

We present a detailed analysis of the evolution of impurity screening cloud in 
the Resonant Level Model (RLM) followed by quenches of system parameters. The 
considered quenches are either local or found to be effectively local. The screening 
cloud is characterized by impurity-bath correlators and by the entanglement of a 
block centered around the impurity with the rest of the system. We consider several 
local quench protocols, such as, quenches from initially decoupled to coupled impurity, 
quenches between different finite couplings of the impurity, and ‘detuning quenches’ 
involving changing the onsite impurity potential away from or to the chemical potential 
of the bath. The relevant correlators and the entanglement display versions of the 
‘light-cone’ effect, since the information about the quench travels through the bath at 
finite speed. At long times (‘inside’ the light cone), the impurity screening cloud relaxes 
exponentially to the final equilibrium structure, with the relaxation rate given by the 
emergent energy scale of impurity screening. Also, snaphots of the time-evolving spatial 
profile of impurity-bath correlators show exponential dependences on space, with the 
length scale given by the screening length. 
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1. Introduction 


The study of real-time evolution of many-body quantum systems driven out of 
equilibrium is a subject of intense current interest, partly motivated by remarkable 
experimental advances in the control of non-equilibrium conditions, both in cold atomic 
systems and solid-state setups. As a result of this interest, many systems and phenomena 
in condensed matter physics are being re-examined in non-equilibrium situations. 
Quantum impurity problems have long been a central paradigm in the study of correlated 
quantum matter. Impurity models have contributed to influential physical ideas, such 
as the Kondo model to the development of the renormalization group and the Caldeira- 
Legget model to the study of decoherence. Quantum impurity models continue today 
to be a source of new physical insights, and appear in a variety of contexts. 

It is therefore not surprising that the study of impurity physics out of equilibrium 
is attracting an increasing amount of attention. Much work has by now been done 
on transport through impurities. More relevant to the present work is the interest 
in real-time dynamics induced by sudden changes (quantum quenches ) of impurity 
Hamiltonians [1-18 . 


Generally, an impurity model involves a zero-dimensional object, often a single site 
or spin, coupled to a bath. The spatial structure of the bath is ignored in much of 
the literature, but this is exactly the aspect that will concern us in this work: namely 
the appearance of a spatial signature of the impurity in the bath. This is the so-called 
‘impurity screening cloud’. Single-impurity models often possess an emergent energy 
scale. The most famous is the celebrated Kondo temperature for the single-impurity 
Kondo model 19,20 , but similar energy scales appear in the single-impurity Anderson 
model 20 221 and the interacting resonant level model |23, 241. There is a length 
scale £ corresponding to this energy scale, which suggests that the bath surrounding 
the impurity is affected differently at distances less than £ from the impurity than at 
larger distances x > £, i.e., there should be a ‘screening cloud’ of radius £ surrounding 
the impurity. Although the impurity screening cloud is difficult to observe directly 
experimentally, calculations have shown that this impurity lcngthscale does in fact 
appear in real-space properties of the bath. The properties (e.g., persistent current 
or conductivity) of a mesoscopic device containing a Kondo or Anderson impurity have 
been found to behave differently if the device size is larger or smaller than the size of 
the Kondo cloud 25-29 . Numerical and variational calculations have found real-space 


properties (e.g., impurity-bath correlation functions, distortion of local density of states, 
entanglement properties, etc) to be different for x < £ and x > £ 30-40 , for Anderson 
and Kondo models and for spin-chain versions of the Kondo model. These recent results 


extend earlier perturbative calculations of real-space structure 41-43 


Perhaps the simplest example of an impurity generating a spatially extended 
screening cloud is that of the resonant level model (RLM), which is a model of non¬ 
interacting spinless fermions. A single ‘impurity’ site (or level) is weakly tunnel-coupled 
to a fermionic bath. The impurity coupling J generates a small energy scale T oc J 2 , 
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and a corresponding large length scale £ oc J 2 . The structure of the equilibrium 
screening cloud in this model was examined in detail in Ref. 44 . The screening cloud 


was characterized primarily through the two-point correlation function between the 
impurity site and bath positions at distance x from the impurity. This impurity-bath 
correlation was found to have different spatial dependence within (x < £) and outside 
(x > £) the cloud. The screening cloud was also described using the entanglement of a 
spatial block containing the impurity with the rest of the bath; the dependence of this 
entanglement on the distance of the block boundary from the impurity is different when 
this distance is smaller than or larger than £. 

In this work, we examine the non-equilibrium response of the screening cloud of 
the RLM after a quantum quench. We start with quenches of the impurity tunnel 
coupling J (section [3]) . This being a local quench, the effect of this change propagates 
through the bath at a finite speed, and after the ‘wavefront’ has passed through a 
certain point, observables related to that position relax to the value corresponding to 
the ground state of the final Hamiltonian. We show that the relaxation of the impurity- 
bath correlator happens exponentially, with the rate of the exponential decay given by 
the final impurity screening energy scale T/ oc -Jj. At any instant, the deviation of the 
correlators from their final value is found to increase exponentially with distance from 
the impurity, with the spatial length scale £. The space- and time-dependence of the 
correlators happens largely, but not completely, through the combination T ft ~ x/f- 
For quenches starting from the impurity-decoupled case, we can explain many features 
of the time evolution of impurity-bath correlators from an analytic calculation in the 
wide-bandwidth approximation. The evolution of the block entanglement also shows 
numerical signatures of the screening energy scale T f , in a somewhat more complicated 
manner. Quenches of the impurity energy (detuning from the bath chemical potential, 
section [4]) results in similar overall behavior as long as the detuning is moderate. If 
the final detuning is nonzero, this provides an additional oscillation frequency in the 
dynamics. Finally, we consider a class of global quenches of bath parameters (section 
[5]), and show how these are equivalent to local quenches for the purposes of screening 
cloud relaxation. 

In the concluding section we provide some discussion and context, and contrast 
several recent studies (lf 3] which are related in spirit to this work. 


2. Model, tight-binding realization, and equilibrium properties 

We consider the resonant level model (RLM), given by the Hamiltonian 

H = £ k cjc k + e d d ] d - J(d t c x=0 + 4 =0 d) . (1) 

k 

Here c*,, c' k (c x , cj.) are the bath fermion operators at momentum k (at position x) and 
d, df are the fermion operators at the impurity site, £k is the dispersion of the bath 
fermions, and J is the hopping strength between the impurity and position x = 0 of 
the bath. The on-site detuning potential £ d is generally tuned to the bath chemical 
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potential, except in Section [4j Here x represents the distance from the impurity. The 
impurity coupling J is assumed to be much smaller than the bandwidth of the dispersion 
£k ; this is the regime where the physics is relatively insensitive to details of the bath 
dispersion. 

At equilibrium, the RLM for small J displays a screening cloud, as examined in 
detail in Ref. [44]. The details depend weakly on the nature of the bath; we will 
concentrate on a one-dimensional bath. Specifically, we will use a closed tight-binding 
chain with 2L — 1 sites, labeled as x = — (L — 1),..., —1, 0,1,..., (L — 1), as shown in 
Figure [jja). We will mostly work with nearest-neighbor (n.n.) hopping in the bath, so 
that the Hamiltonian is 

L—l 

H = - ^2 ( c I c m + c l+i c i) + £ dd ] d - J (cjd + d ] c 0 ^ . (2) 

i=-L -\-1 


The site label L is identified with site — L + 1 as usual for periodic boundary conditions. 
We take the bath n.n. hopping to be unity, so that energies (time) is measured in units 
of the bath n.n. hopping strength (inverse of the bath n.n. hopping strength). The 
dispersion is then e& = —2 cos fc, with k = 27rn/(2L — 1) and n = — L + 1, ...,L — 1. 
In Section [5j we also modify the bath by adding next-nearest-neighbor hoppings to the 
tight-binding chain. We will also mostly consider half-filling, i.e., L particles, since the 
total system including the impurity has 2 L sites. The bandwidth in these units is 4, we 
are thus interested in the parameter regime J <C 4. In practice, we will restrict the final 
J values to be < 0.5. 

The equilibrium screening cloud around the impurity is well described by the single 
particle impurity-bath correlator (d)c x ) 44 . The zero-temperature equilibrium spatial 


structure of ( d)c x ) features a rapid oscillatory behaviour with wavelength given by the 
Fermi momentum k F , convoluted with an envelope function that crosses over from a 
oc In x behaviour at small x, to an algebraic decay oc x _1 at large x. The length scale 
£ at which this crossover happens, i.e., the ‘size’ of the impurity screening cloud, is 


x 


£ = Vf/ttJ 2 p(£f), scaling as J~ 2 with the impurity-bath coupling. Here v F the Fermi 
velocity and p(e F ) the normalized single-particle density of states [J dup(u) = 1] of the 
bath at the Fermi energy e F . With our units, at half-filling vp = 2 and p(e F ) = 1/(27t), 
so that £ = 4/J 2 . 

In Ref. (44), the equilibrium screening cloud was also characterized through a study 
of the entanglement entropy, motivated by similar characterizations in Refs. [39||40| > |45] . 
For this purpose, one considers the entanglement entropy Sj^ = —Ti 'A t [pA t In pA t ] of the 
subsystem A 1 of total size 2 i centered around (and containing) the impurity. Here pA t 
is the reduced density matrix of the subsystem An obtained by tracing out the degrees 
of freedom lying outside the subsystem, starting from the wavefunction (in this case the 
ground state wavefunction) of the full system. The so-called impurity entropy, defined 


as 


cM _ cM _ cM 

°imp ~ °J °J=0i 


( 3 ) 
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(a) 





Figure 1. (a) Geometry of tight-binding model ([2|. (b) Typical behavior of real and 

imaginary parts of (d'c x ) at x = 24, after a quench from J = 0.8 to J = 0.3. The two 
dashed vertical lines and the associated cartoons indicate the time for the wavefront 
to reach the point x (~ x/2), and the time for finite-size effects (~ 2L+ f~ x ). (c) The 
propagating wavefront, seen visually by comparing the time evolution of ( d'c x ) for 
different x. 


was shown to have, as a function of the subsystem size £, a crossover to a an algebraic 
oc £ _1 decay at large £ around the same characteristic length scale £ 

In the following we analyse the non-equilibrium dynamics of the screening cloud 
after a quantum quench by monitoring the evolution of ( dfc x ) and as a function 
of time. For a system initialised in the ground-state of the Hamiltonian for t < 0, we 
consider local quenches, where either J or are changed abruptly at t — 0 to their final 
values, and also a quench in the global parameters corresponding to an abrupt change 
of the dispersion relation of the bath. For a bath with an infinite number of degrees 
of freedom, local correlators are expected to relax to the local equilibrium values for 
asymptotic long times after a local quench. For these cases we analyse the way the new 
equilibrium form of the screening cloud is reached. 


44 


3. Quenches of the impurity coupling J 


In this section we will consider quenches of the impurity-bath coupling J. The system 
starts in the ground state for J = Ji and, after the J* —* Jf quench, evolves under the 
final Hamiltonian with J = Jf. In subsections 3.1, 3.2 3.3 we consider the evolution of 
the equal-time impurity-bath correlation function 

C(x,t) = (d\t)c x {t)) . (4) 


In subsection 3.4 we consider the evolution of the block entanglement entropy. 
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3.1. Time evolution of impurity-bath correlators: general features 

In figures [ljb) and [l|c) we display some general features of the time evolution of the 
impurity-bath correlator, highlighting effects of locality, information propagation, and 
the appearance of infinite-size quench features in finite-size systems. 

After a local quench, information about the change in the system should propagate 
out from the point of disturbance with speed vp. This effect has been discussed widely in 
recent times in terms of the so-called Lieb-Robinson bound and the ‘light-cone’ effect, 
for both local and global quenches [l]-3][46]-[58|. For local quenches, one typical way 
of formulating this causality effect is that a local observable at position x will be 
substantially affected by a local quench at position y only in the space-time region 
t > \x — y\/v (‘inside’ the light cone). Here v is the characteristic speed. 

The equal-time impurity-bath correlator C(x,t) = ( d'c x ) is not a local observable. 

Thus, it evolves in time also outside the light cone, t < x/vf , but this part of the 
dynamics is entirely due to the d operator and hence relatively simple. The dynamics 
is dramatically affected at t ~ x/vp, when the ‘wavefront’ reaches the point x. This is 
shown in the examples of Figure [jjb), in particular the top left schematic. The panels 
in Figure [IJc) visually display the increasing time required for the wavefront or signal 
to reach farther points. 

Within the light-cone, i.e., for t > x/vp, the correlator decays approximately 
(upto ~ L^ 1 effects) toward its final ground-state value, C^ q (x,t). The reasoning 
behind this is that a local quench is a L _1 effect and therefore does not inject a 
thermodynamically significant energy into the system, hence a large enough system 
should relax approximately to its ground state. Figure [l](b) also shows that finite-size 
effects are strongly felt when the signal or wavefront traverses the other end of the system 
and then reaches point x (top right schematic). Until this point of time, t & (2 L — x)/vf, 
the correlator (d'c x ) is approximately L-independent, thus by analyzing this region we 
essentially access the behavior of (dfc x ) in an infinite system. By increasing the ring size 
L, we get larger windows of time where the infinite-system dynamics can be accessed. 

The details of the dynamics depend sensitively on the parameters, filling, and the 
site x. Even for the half-filling case (where the behavior is cleanest, also in equilibrium), 
the real and imaginary parts of C(x,t) show different dynamics for odd and even x. 

For all x, the imaginary part starts from zero and (for large enough systems) goes to 
zero at long times. For even distances x , e.g., Figure [jjb), the imaginary part shows 
strong finite-size effects and for large enough systems is almost time-independent. For 
odd x (not shown), the behavior is reversed: Re[C(x,£)] shows strong L-dependence 
while ImjC^x, £)] does not. In order to avoid these system- and ^-dependent details, we 
will mostly show absolute values, \C(x,t)\ or | C(x,t) — C^ q (x,t)\. 

A local quantity such as (c\.c x ) or (c\.c x+ f) shows a clearer light-cone effect, i.e., 
stays constant for t < x/vp until the wavefront reaches x. Such data is not shown in 
this work, as our focus is the impurity screening cloud, which is better embodied in the 
impurity-bath correlator C(x,t ) = {dfc x ). Also, Ref. [l] has discussed the appearance 
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t — x 


Figure 2. Quenches from an initially decoupled impurity (Ji = 0). The solid lines 
are tight-binding numerical results for sizes L ~ 1000. Dashed lines are the infinite- 
bandwidth infinite-size results of Eq. ([7| . (a) Snapshots of the rescaled impurity-bath 
correlation function C(x,t) as a function of x = x/£f for different times. f(x) is the 
scaling function obtained in Ref. [44]. (b,c) Time evolution of the difference between 
C(x,t) and its equilibrium value C^ q (x), at several positions x. We plot against 
the rescaled and shifted time variable t — x. (d) Snapshots of the spatial profile of 


C(x, t ) — C^ q (x) 


at time t = 100. 


or absence of the light-cone effect in different unequal -time correlators. 


3.2. Starting from decoupled impurity: Ji = 0 


We first consider an initially decoupled impurity level which is suddenly coupled at time 
t = 0 by changing the impurity hopping parameter J from Jj = 0 to finite Jf. The 
initial correlation functions are thus given by 

(4 c fc')o = dkk'Q{k F - \k\)] (dfd )o = n d \ (d ] c x ) = 0 (5) 


where n d is the initial occupancy of the impurity level. In the final equilibrium state, 
with J = Jf 0, the impurity is screened over a finite distance £/ = vf/Lj, with 
Tf = 7Tjjp(£ F ). 

In these conditions the evolution of the equal-time two-point correlator (df (t)ci(t)) 
is amenable to an analytic treatment along the same lines of Ref. jT]. The approach is 
detailed in the Appendix. We will express results in terms of the rescaled correlator 


J 

7 T 2 V F 


C(x, t ). 


C(x, t) 


( 6 ) 
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In the thermodynamic limit and the large-bandwidth limit, we obtain the following 
expression for the time evolution of the rescaled impurity-bath correlator: 

„ 0 1 — e ~i(e-i)i 


C(x, t) = it / de- 


cos, 


X ( £ + k F 


(e-i) 


6(t — x) [ine ikFX ] 


dee 1 


/ £ 2 + 1 
' —oo ° 1 


3 — i(e—i)t. 


0 i(e+i)(t-x) 


( 7 ) 


Here x — x/£f, t — T ft, and kF = ffkF are dimensionless versions of distance, time and 
Fermi momentum, obtained by scaling with the length scale ff and energy scale Tf of 
the impurity screening cloud corresponding to the final impurity coupling Jf. For our 
lattice Hamiltonian ([2]) at half-filling, kp = tt/2. 

We comment on several features in this expression: 


• There is a clear and natural separation between behaviors for t < x and for t > x, 
i.e., inside and outside the light cone. The first term shows that the evolution is not 
frozen for t < x (outside the light-cone). As discussed in the previous section, this 
is expected in the impurity-bath correlator because of the evolution of the impurity 
operator df{t) 

• The impurity screening cloud (associated with the ground state of the final 
Hamiltonian) sets time- and length- scales for some of the dynamics, since time 
and position appear in the combinations t — F ft and x = x/ff in the expression. 
In particular, the non-oscillatory part of the time-dependence factorizes out of the 
integrals and is of the form e _t , which shows that the long-time decay to the final 
value is exponential in time with time constant 1 /F f. 

• The time-dependent correlator displays kp oscillations in space, 
e ik F x _ e ik F x, j US £ like the equilibrium behavior of the correlator 

• In the t > x part (inside the light-cone), a non-oscillatory exponential spatial 
dependence e x can be factored out of the integral. Thus, the deviation from the 
final equilibrium value increases exponentially with distance, in the spatial region 
in which the wavefront has already propagated through. 


ue to factors like 
44 . 


The long-time limit (t —» oo) of this expression reduces to the equilibrium expression 
for the impurity-bath correlator derived in Ref. 44 , showing that for large sizes 
the impurity-bath correlator reduces to the final equilibrium profile. 


• There are no temporal oscillations with frequency = 2 (twice the bath hopping 
or half the bath bandwidth). In calculations with the tight-binding model [2j 
oscillations at this frequency appear, as in Figure [T|(b,c) ; this represents the energy 
width of occupied single-particle states. However, since the expression is obtained 
in the infinite-bandwidth limit, the bandwidth or hopping parameter of the bath 
do not appear as an oscillation time-scale. 

Figure [ 2 ] shows the impurity-bath correlation function C(x,t ) after quenches 
starting from the decoupled case, Jj = 0, comparing numerical tight-binding results 
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using Hamiltonian ([2]) against predictions from Eq. ([7]) above derived in the large-size 
large-bandwidth limit. 

Figure a) shows snapshots of the spatial prohle of \C(x,t)\, calculated for the 
tight-binding Hamiltonian at different time instants t after the quench. The spatial 
prohle is seen to approach the final equilibrium curve at large times. The (final) 
equilibrium prohle approximates the analytic equilibrium curve f(x), but does not match 
it exactly for a single finite system size and single Jf value, as discussed in detail in 


Ref. 44 


To highlight the approach to equilibrium, in Figures |2](b,c,d) we show the correlator 
with the final equilibrium value subtracted oh, i.e, C(x,t ) — C^ q (x). Motivated by the 
light-cone structure, we plot this as a function of the variable t — x. The behavior 
outside the light-cone ( t < x) and inside the light-cone ( t > x) look sharply different, in 
this quantity. The quantity is plotted for different bath positions x (Figure [2](b) ), and 
for different final coupling Jf at the same bath position x = 100 (Figure [2](c)). In the 
log-linear plots, the decay for t > x appear as straight lines with the same slope in each 
of these cases, showing that the decay to the final value inside the light cone happens 
exponentially with the time scale associated with the (final) impurity screening cloud. 

Figure [2](d) plots spatial snapshots of the correlator at the same time instant, 
t = 100. We choose to plot it against the same variable t — x used for showing the 
time evolution in Figures [2j)b,c) . The similarity of paneled) with the plots in i>. c) 
highlights the fact that the dynamical features are dominanted by the quantity t — x. 
In Figure |2](d) we clearly see an exponential increase of the deviation from the final 
value, with increasing distance. This is consistent with the e x factor in the infinite- 
bandwidth expression, discussed above. There are small e lkpX oscillations, but they are 
negligible everywhere except at small values of x. These oscillations appear in both the 
tight-binding data and in the infinite-bandwidth equation, and show strong finite-size 
effects. 

In each panel of Figure [2j we show tight-binding results (solid curves) together 
with the prediction Eq. [ 7 ] (dashed curves). For the regime under discussion, the infinite- 
bandwidth continuum approximation captures well the essential features of the impurity 
cloud dynamics as seen through the C(x,t) correlator. The agreement is generally 
excellent. Points of disagreement are the temporal oscillations of frequency 2 (half the 
bandwidth) appearing in the lattice data but not in the infinite-bandwidth analytics, 
and the difference between the e lkFX oscillations, which are anyway very small. 


3.3. Quench from finite Ji to finite Jf 

We now consider quenches from finite Ji ^ 0 to a different finite value Jf. The initial 
correlators now no longer have a simple form like Eq. (|5]) for the Ji = 0 case, but rather 
they describe a screening cloud with characteristic length Q oc J* -2 . 

After the quench, the correlator C(x,t) relaxes in an analogous way to the new 
equilibrium prohle corresponding to the length ff oc Jf~ 2 ■ As the light cone spreads, 
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Figure 3. Quenches from finite Ji to finite J/. (a) Relaxation of the impurity- 

bath correlation function, for x = 50. In (b,c), the results are averaged over small 
time intervals At to average out the strong oscillations and emphasize the exponential 
decay with rate Tf. In panels (b) and (c), At = 20 and At = 12 respectively. 


C(x,t) within the light cone (x < v F t) at time t will relax exponentially to the final 
profile. 

In principle, an integral expression for the correlation function C(x,t ) could be 
derived analytically using a similar method to the one given in the appendix for J, : = 0. 
However, such an expression would involve additional integrals, because the initial 
C(x,0) is already nontrivial 44 . We therefore restrict to showing numerical results 


from the tight-binding model. For Ji ^ 0, the temporal oscillations of C(x,t) (with 
frequency e F ) are much stronger than in the J t = 0 case, as seen in Figure[3|a). To focus 
on the exponential relaxation, in Figures |3](b) and[3|c) we average out the oscillations 
by plotting the running average of C(x,t), i.e., each point is an average of C(x,t ) taken 
over a small time window At around time t. For a variety of J, and Jf values, the 
relaxation is seen to be always exponential, with the time scale corresponding to the 
final impurity screening scale Vf = nJjp(£ F ). 


3 . 4 . Entanglement entropy 

We now consider the time evolution of the entanglement entropy (t) between a block 
A containing the impurity at the center, and the rest of the system, which we call B. 
The geometry is shown in Figure |4](a). The block A contains a total of 2£ sites including 
the impurity, so that the boundary between the blocks is at distance £ from the impurity. 

The entanglement entropy in a free-fermion system can be obtained from the matrix 
of all two-point correlators within one of the blocks |44[[59||60|: if zy’s are the eigenvalues 
of this matrix, then S = — yT[zyliizy + (1 — z/j)ln(l — uf)]. In the present case, we 
calculate the entanglement at instant t from the matrix of equal-time two-point 

correlators at time t. 


Ref. 59 has studied the entanglement evolution after a particular local quench, for 


a block containing the quench position at the center. (The strength of a single bond of 
the chain is quenched.) The results are intriguing and, to the best of our knowledge, 
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0 200 400 

Time t 


100 1000 

Time t 


Figure 4. (a) Geometry for defining block entanglement between A and B partitions, 
(b) Time evolution of entanglement entropy, with l = 40, after a quench from Jj = 0 
to two different final Jf. The dynamics starts at t = 1/vf = 20, after the signal 
reaches bourndary between A and B blocks. The insets zoom around this instant. 
The dynamics clearly is faster for larger Jf. The entanglement entropy approaches the 
final ground-state value (dashed horizontal lines) at long times, (c) Effect of Jf on the 
relaxation explored in detail, for quenches starting from Ji = 0 and from Ji = 0.5. The 
left panels (log-linear scale) highlight the initial ~e _r/t behavior. The right panels 
(same data, log-log scale) show that the long-time behavior is an approximate power 
law, with drifting exponent that suggests a logarithmic correction as in Ref. 


59 


not yet completely understood: the block entanglement entropy decays to the final 
equilibrium value with a power law modified logarithmically. This is accompanied by 
In t behaviors of some of the eigenvalues ig of the correlation matrix. Our situation is 
broadly similar, but we are particularly interested in the role of the impurity screening 
scale in the dynamics. 

Figure |d](b) shows the behavior of S^(t) after a quench starting at Jj = 0. (The 
overall behavior for other Ji values is very similar.) The entanglement between the 
blocks stays frozen at the initial value until the information wavefront reaches the 
boundary between the blocks, i.e., until t ~ 1/vf- This may seem surprising because 
many correlators in the correlation matrix of A have dynamics at smaller times (as 
seen in previous sections), but it is physically reasonable because 5^(t) measures the 
entanglement between A and B blocks, and the B block cannot know about the local 
quench until this time. After a sharp reaction shortly after t > £/vp, the entanglement 
relaxes toward the final ground state value, shown in Figure |4](b) with dashed horizontal 
lines. The overall form of (t) is broadly similar to that seen in the local quench of 
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t > 0 





t t 


Figure 5. (a) Schematic showing a quench from finite to zero detuning. (b,c) 

Relaxation of impurity-bath correlators for quenches to zero detuning. (d,e) Evolution 
of impurity-bath correlators for quenches to finite detuning. In order to compare 
frequencies, in these panels we plot against time without rescaling. 


Ref. 59 . The response and relaxation behavior are clearly faster at larger Jf, consistent 
with the fact that the energy scale of the screening cloud (T f) increases with Jf. 

The relaxation behavior is examined in more detail in Figure [4j^c) . There are two 
time regimes, separated by a time scale proportional to T/ Jf 2 . (For Jf = 0.1, only 
the first regime is visible.) The initial relaxation regime is exponential, with the time 
constant given by T^ 1 . This feature is related to the physics of impurity screening and 
does not appear in the quench of Ref. |59|. At later times, the relaxation is slower; 
the data for accessible sizes is consistent with a power law but the best-fit exponents 
vary. This suggests that the exact behavior in this long-time regime may well have 
logarithmic corrections as in Ref. 59 . The time evolution of the eigenvalues zy also 


shows logarithmic evolution of some eigenvalues as in Ref. 59 ; a full understanding is 
lacking at present, but presumably the In t features are not related to impurity screening 
physics. 


4. Quenches of the impurity detuning e^ 

In this section we study the dynamics after quenching the detuning, i.e., the onsite 
energy bias of the impurity site, denoted by in Equations (JTJ) and (j2]). The detuning 
energy is measured with respect to the Fermi energy, as shown in Figure [5|a). We 
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consider both quenches from an initially detuned impurity potential to the zero detuning 
case, £ d 0 to £ d = 0, and vice versa. The first situation is depicted in Figure |H](a). 
The impurity coupling J is kept unchanged in this quench protocol. 

For the ground state, tuning of the on-site impurity potential changes the structure 
of the impurity cloud at small distances (within the screening cloud, x < £) and 
eventually destroys the logarithmic dependence oc In x of the correlation function (dfei) 
in this spatial region 44 . However, the effect is significant only at large £d [44]. Here, 


we consider moderate values of £d, of the order of the bath bandwidth. 

Figures |H](b,c) show quenches of the detuning to a final zero value, £ d = 0. The 
relaxation of correlators to the final ground state value is exponential, with the relaxation 
rate set by the screening cloud scale, T/ = J 2 / 2, as before. The data shows the usual 
half-bandwidth frequency oscillations. In Figure [5|c) the oscillation periods appear to 
be differentfor different J, but this is only due to the J-dependent scaling of the time 


axis. 


In Figures |5[d,e) we show the dynamics after quenches to 0(1) values of £ d . The 
decay to the final values occur exponentially, with the same time scale ITT, consistent 


44 


However, 


with the fact that moderate £ d values do not destroy the screening cloud 
there are now strong oscillations with freuquency £ d . To highlight the frequencies, we 
use axes without dividing by the screeing scale. The correlators now show oscillations 
with two frequencies — one equal to half the bandwidth and one equal to £ d . In Figure 
[H](d) the new frequency £ d — 1 (prominient oscillations at smaller times) is half that of 
the smaller oscillations visible at later times. In |5](e) the oscillation frequency visibly 
varies with £ d . 

At extremely large final detuning values, £ d 1 (not shown), the relaxation appears 
to be power-law-like rather than exponential, for accessible system sizes. Since this value 
is larger than the bandwidth, this is beyond the regime of usual impurity phsyics, and 
is consistent with the absence of impurity screening. 


5. Global quenches with local effects 

In this section we consider a global quench which behaves effectively like a local quench, 
in the sense that the correlators relax to the ground state impurity screening cloud 
corresponding to the final Hamiltonian. We consider a quench that changes the 
dispersion of the bath, through a next-nearest neighbor hopping with opposite sign 
compared to the nearest-neighbor hopping. The Hamiltonian is now 

L—l L—l 

H = - (4^+1 + 4+iq) + j nn ( c ^ Ci + 2 + c ^+2 c *) 

Z—— L~\~ 1 Z—— L -|-1 

+ £dd)d — J + c^coj . (8) 

This corresponds to a dispersion relation £k = —2 J cos (k) + 2J^n cos(2 k), as shown in 
the top cartoon in Figure |6](a). 
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t — x 


t — x 


Figure 6. (a) Schematic showing a quench from finite Jnn to Jnn = 0. (b,c) 

f 

Relaxation of impurity-bath correlators, for various quenches ending at J J NN = 0. 
Panel (b) uses various initial Jnn values and panel (c) uses various J values. The 
behavior is remarkably similar to the case of detuning quenches. 


The equilibrium (zero temperature) effects of such a modification to the bath 
dispersion was explored in some detail in Ref. 44 . At any filling fraction, there are 


three regimes of behavior depending on the value of Jnn- For half-filling these behaviors 
are as follows: For Jnn > 1, the bath has two pairs of Fermi points rather than a single 
pair; the impurity screening cloud has the same form except that the density of states 
now has contributions from four Fermi points. The Jnn — 1 case is special as there 
is single Fermi ‘touching’ point instead of a pair of Fermi crossing points. In this case 
the correlators show a non-Fermi-liquid character for distances exceeding the screening 
length. Finally, the 0 < Jnn < 1 regime is very similar to the case of the un-modified 
bath, with two Fermi points. 

In the absence of the impurity, the single-particle Hamiltonian with and without 
the Jnn term has the same eigenstates for each momentum, only the eigen-energies 
are affected. For any Jnn £ [0,1], the same single-particle momenta k e [—|, |] 
constitute the ground-state Fermi sea at half filling, so that an abrupt change of J N n 
in the absence of an impurity would leave the system stationary and not have any non¬ 
equilibrium effects. Such a quench has an effect only due to the presence of an impurity, 
so that it is effectively a local quench. For present purposes we consider half-filling and 
0 < Jnn < 1 - 

Figure [6] shows data for quenches from nonzero Jnn to Jnn = 0, with the detuning 
parameter set to zero: £d = 0. The initial equilibrium state has two Fermi points as in 
the case of an un-modified bath; however the chemical potential is no longer at zero for 
half-filling, so that the detuning compared to the Fermi surface is initially nonzero. Such 
a quench is thus effectively a detuning quench of the impurity. In fact, the results for 
the time evolution of impurity-bath correlators shown in Figure [6](b,c) are remarkably 
similar to those in Figures[5](b,c), where the detuning was quenched from nonzero values 
to zero. The change in the dispersion modifies vp and p(sf)', this also has an effect on 
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the dynamics of impurity-bath correlators, but the dominant effect is clearly the change 
of detuning between Fermi level and impurity energy. 

We have also performed quenches from J^n = 0 to Jnn > 0 (not shown). 
Repeating the argument above, these are similar to detuning quenches from zero to 
nonzero detuning. Accordingly, the time evolution of C(x,t) = (d'c x ) is very similar to 
the detuning quench results of Figures [5](d,e). 
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6. Discussion 


We have presented a study of the dynamics of the impurity screening cloud in the 
resonant level model following local (or effectively local) quenches. We have focused 
mostly on the time-dependence of the equal-time impurity-bath correlators C(x,t ) = 


(d' (t)c x (t)), which characterizes the screening cloud at equilibrium 44 


Due to the finite speed of propagation in the bath, the effect of the quench is felt 
at point x only at times t > x/vp■ In terms of the correlator C(x,t), this means that 
for t < x/vf the dynamics is due only to the evolution of the d) operator. Thus, C(x, t ) 
has relatively simple evolution for t < x/vp. After the information wavefront has passed 
point x, the correlator relaxes to the ground state of the final Hamiltonian, provided 
the system is large enough. We have outlined how the behavior is modified for finite- 
size baths, and how the infinite-system quench properties can be extracted from time 
evolution data on finite-size systems. 

The relaxation to the final ground state value is found to happen exponentially , 
with the rate given by the energy scale T f associated with impurity screening in the 


final Hamiltonian. This is the scale of broadening of the impurity spectral function 44 
analogous to the Kondo temperature scale in the single-impurity Hondo model. 


We 

axation with this scale after the wavefront passes the point at 
C(x, t ) — Cf q (x) 


have highlighted the re 
time x/vp, by plotting C(x, t) — C^ q {x) against the rescaled and shifted time variable 
Tf(t — — t — x for a variety of quenches, with T f determined by the final impurity 

coupling ,Jf. The relaxation curves plotted in this way look broadly similar for a range 
of quenches. For the case of a quench starting from Jj = 0, we derived and presented an 
explicit expression for C(x,t), which analytically expresses these physical effects, e.g., 
through explicit appearance of &(t — x) factors signaling the information wavefront. 

We have also described the effect of a detuning quench on the equal-time correlators. 
For moderate detunings, the exponential relaxation persists, and if there is a final 
nonzero detuning with respect to the chemical potential, this provides an additional 
oscillation frequency. A quench of the bath involving nearest-neighbor hoppings reduces 
effectively to a local quech as long as the number of Fermi points is preserved. 

The evolution of the entanglement entropy of a block centered around the impurity 
shows a clear light-cone effect (no dynamics before t ~ x/vf). For t > x/vf, we find 
two relaxational regimes. The first regime shows exponential relaxation with the rate 
T f. The longer-time regime is consistent with a power-law relaxation with logarithmic 
corrections, as in the case of Ref. 


59 


The dynamics of impurity screening clouds has been the subject of a few recent 
studies [l]-j3], complementary to ours. We briefly outline relevant parts of these results 
below and compare with our results. 

In Ref. [l], the authors have studied the Kondo model at the so-called Toulouse 
point 61 . At the Toulouse point, the Kondo model reduces to the non-interacting 


resonant level model, and hence various analytical results can be derived in the case 
when one starts from zero coupling. The authors focus on spin-spin correlators of the 
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Kondo model, which translate to different fermionic correlators from the ones occurring 
naturally in this work. Focusing on unequal-time correlators, Ref. [l] distinguishes 
correlators which show a strict light cone from those which do not. 

In Ref. |2|, quenches in the isotropic Kondo model are studied using time-dependent 
NRG and perturbation theory, for both ferromagnetic and antiferromagnetic Kondo 
couplings. The time evolution of equal-time impurity-bath spin-spin correlators are 
studied. As in our case, the correlators considered do not have a strict light cone. 

Ref. (3] addresses several issues similar to the ones addressed in the present work. 
In this work, the dynamics of the Anderson model is studied with a focus on the Kondo 
regime. Two types of equal-time impurity-bath correlators are studied: spin-spin and 
density-density correlators. The behavior for t < xjvp and t > x/vp are characterized 
in some detail, and the exponential decay for t > x/vp is discussed. The decay rate 
appears to be more complicated than in our case and weakly dependent on position; it 
is not completely clear whether this is due to system size issues. 

The present study has focused on the noninteracting RLM and the dynamics of 
correlators naturally suitable for describing the RLM screening cloud. The absence 
of interactions has allowed us to numerically treat relatively large systems, making 
possible an unambiguous identification of the time scale of relaxation. This aspect is 
more difficult with sizes accessible for interacting systems |3j. In addition, our analytical 
result for the Jj = 0 case explicitly shows ~e _r A relaxation of the equal-time impurity- 
bath correlators. Based on the present results, it is natural to conjecture that the 
relaxation rate is given by the characteristic equilibrium impurity energy scale for any 
impurity model characterized by an emergent energy scale. For example, the Kondo 
temperature scale should provide the relaxation rate for the antiferromagnetic Kondo 
case of Ref. |2|, but the details of relaxation are not yet clarified in the literature. 

The previous studies concentrate on initially decoupled impurities, i.e., cases where 
there is no nontrivial screening cloud in the initial state. In section 33 we have presented 
numerical results for finite initial coupling J*. Even though the relaxation is governed 
solely by the final impurity energy scale, the structure of the initial screening cloud 
induces strong oscillations superposed on the relaxation. 

Local quenches are being currently studied from a variety of perspectives. The 
present work, together with the related literature discussed above, addresses the role of 
the emergent energy scale in the dynamics, in cases where the local feature possesses 
such a scale. Several aspects remain poorly understood and clearly deserve further 
investigation and explanation, such as the behavior of the block entanglement entropy 
and entanglement spectrum, the role of finite bandwidths, etc. 
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7. Details of the derivation 


Eq. (Jt|) , reported in the main text, gives the time evolution of impurity-bath correlators 
in the infinite-bandwidth and infinite-bath-size limit. In this appendix we briefly outline 
the derivation of this expression. 

The evolution of the system after the quench at time t — 0 is governed by the final 
Hamiltonian given in Eq.Q which in its eigen-operator basis a £ can be written in the 
form H = Yl,E £a \ a E■ The new set of operators a e , ai can be related with the original 
Cfc and d by imposing the commutation relations [H, a e ] = — ea £ . This yields 

a e = A ed d + A £ kCk, a\ = A* d df + ^ A* Ek c\ (9) 

k k 

with 


A ed = AI* = 


ed 




A £k = AU = 


2 L +1 r 2 + e 2 ’ 

,]' 1 


ek 


2 vi 


\/2L + 1 £ — £fc V 2T + 1 T 2 + s 2 

Note that the unitarity of the transformation also implies that 


( 10 ) 

( 11 ) 


Ck 


^ ^ A-zk&ei d ^ ^ A £c l& £ . 
£ £ 


( 12 ) 


The time evolution of the operators df{t) and c r (t) can be most easily computed in 
the basis where H is diagonal. This leads to the expressions: 

d\t) = V m Ed | 2 e“*d+(0) + J2 A ^ e<et 4 ( 0 ). (13) 

£ £,k 


C x (t ) = 


V2L + 1 . 


^ ' A £k A £d e 


ikx ~—i£t 


d(0) 


k,£ 


^ ^ A £ kA £ ^c 


ikx „—i£t 


Cfc'(0) 




(14) 


relating the time evolved operators to the operators at t — 0. In order to obtain a useful 
expression for the correlator (d'(t)c x (t)), we now introduce two approximations. The 
first is the thermodynamic limit where an infinite chain is considered. The second is 
the wide-band approximation where we assume that the energy scales of interest are 
much smaller than the bandwidth of the chain. This amounts to considering that the 
density of states of the bath is constant. With these two approximations in place, the 
summations over momentum can be done assuming a linear dispersion relation around 
the two Fermi points: e k = Vp(k — kp) and e k = ~vp(k + kp). Thus for any function 
f(k) of momentum we have: 

2 ^ + y -» J [f( kF + v F l£ k) + f(-k F -vJ 1 e k )] (15) 

However, some of the coefficients of the unitary transformation have denominators of 
the form which have to be treated with care when e ~ e k . By analysing the discrete 
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version of the eigenvalue equation of the H : 

1 

£-£ k 

one arrives at the prescription (see also Ref. [l] and references therein) 


J 2 


£ = 


2L + 1 


E 


->• P 


£~£k 


£ ~ £k 


+ 7r-(5(e-e fc ) 


(16) 


(17) 


where P stands for the principle value. Within the considered approximations, the 
summations in Eqs.(13, 14) can thus be replaced by integrals and partially done by 
contour deformation in the complex plane. An example of such a procedure for a typical 
term arising in the computation of c x (t) is: 


V ^ e,k J 


-,—iet 


d£ 


£ 2 + A 2 


E 

s=±l ' 


d£ 


2nvi 


^ e is{k F +v F 1 e k )x 


X 


-“- J + _ £ k) 

£ ^ k J 

= -i—e^Oit-x^e-^ 
Vf 


(18) 


where we assumed x > 0 for simplicity. The pole structure of the expression leads to 
the theta-function term 9(t — x) in Eq.Q that signals the passage of the wave front 
through the rescaled position x. 
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